Parallel Iterative Solution of Sparse Linear Systems Using Orderings from Graph Coloring Heuristics

The eeciency of a parallel implementation of the conjugate gradient method precon-ditioned by an incomplete Cholesky factorization can vary dramatically depending on the column ordering chosen. One method to minimize the number of major parallel steps is to choose an ordering based on a coloring of the symmetric graph representing the nonzero adjacency structure of the matrix. In this paper, we compare the performance of the preconditioned conjugate gradient method using these coloring orderings with a number of standard orderings on matrices arising from applications in structural engineering. Because optimal colorings for these systems may not be a priori known, we employ several graph coloring heuristics to obtain consistent colorings. Based on lower bounds obtained from the local structure of these systems, we nd that the colorings determined by these heuristics are nearly optimal. For these problems, we nd that the increase in parallelism aaorded by the coloring-based orderings more than oosets any increase in the number of iterations required for the convergence of the conjugate gradient algorithm.